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The branching fractions of the exclusive decays K'^'y and — > K*^'y are measured 

from a sample of (22.74 ± 0.36) x lO** BB decays collected with the BABAR detector at the PEP 
II asymmetric e+e" collider. We find B{B° K'^-y) = (4.23 ± 0.40(stat.) ± 0.22(sys.)) x 10"^ 
B{B+ K'+'y) = (3.83 ± 0.62(stat.) ± 0.22(sys.)) x 10"^ and constrain the CP-violating charge 
asymmetry to be -0.170 < Acp{B K*^) < 0.082 at 90% C.L. 

PACS numbers: 12.15.Hh, ll.30.Er, 13.25.Hw 



In the Standard Model (SM) the exclusive decays 
B — > if *7 proceed dominantly by the electromagnetic 
loop "penguin" transition b sj. Many extensions 
of the SM provide new virtual high-mass fermions and 
bosons that can appear in the loop, causing deviations in 
the inclusive rate for 6 ^ 57 [|l| . The sensitivity of the ex- 
clusive rates to these effects is limited by the uncertainty 
in the SM calculation. However, there has been con- 
siderable progress recently |Q. The precision measure- 
ment of the exclusive branching fractions B{B'^ — > K*'^j), 
B{B^ K*~^'y) is needed to test and improve these cal- 
culations. The non-SM processes can also interfere with 
the SM decay to cause CP- violating charge asymmetries 
at a level as high as 20% ||]. The CP- violating charge 
asymmetry from SM contributions alone is expected to 
be < 1%. 

In this letter, measurements of the exclusive branch- 
ing fractions, B{B° K*°j) in the K*" K+n-, 
K°Tr° modes, and S(B+ K*+j) in the K*+ K+tt°, 
K^'K'^ modes with tt+tt^, are presented. Here 

K* refers to the if* (892) resonance and the charge con- 
jugate decays are implied unless otherwise stated. The 
K*o K+TT- and K*+ K+Tr°, K°Tr+ modes are used 
to search for CP-violating charge asymmetries. 

The data were collected with the BABAR detector |Q] at 
the PEP-II asymmetric e+(3.1 GeV) - e-(9 GeV) storage 
ring 1^. The results in this paper are based upon an in- 
tegrated luminosity of 20.7 fb~^ of data corresponding to 
(22.74±0.36) X 10*^ BB meson pairs recorded at the T{AS) 
resonance ( "on-resonance" ) and 2.6 fb^^ at 40 MeV below 
this energy ("off-resonance"). 

We use Monte Carlo simulations of the BABAR detec- 
tor based on GEANT 3.21 Q to optimize our selection 
criteria and to determine signal efficiencies. Events taken 
from random triggers are used to measure the beam back- 
grounds. These simulations take account of varying de- 
tector conditions and beam backgrounds. 

The selection criteria for this analysis are optimized 
to maximize S'^/{S + B) where S is the number of sig- 
nal candidates expected, assuming the central values 
of the previous measurement B{B'^ K*^^^ B~^ — > 
K*+-/) = {A.55t"all{stat.)±0Misys.),3.76toilistat.)± 
0.28{sys.)) X 10^^ Q, and B is the expected number 
of background candidates determined from Monte Carlo 
and confirmed with off-resonance data. Quantities are 
computed in both the laboratory frame and the center- 



of-mass frame of the e^e^ system. Those computed in 
the center-of-mass frame are denoted by an asterisk; e.g. 
^beam — 5.29 GcV is the on-resonance energy of the e+ 
and e~ beams. 

We require a high-energy radiative photon candidate 
with energy 1.5 < < 4.5 GcV in the laboratory frame 
and 2.30 < E* < 2.85 GeV in the center-of-mass frame. 
A photon candidate is defined as a localized energy max- 
imum Q in the calorimeter acceptance —0.74 < cos6' < 
0.93, where 6 is the polar angle to the detector axis. It 
must be isolated by 25 cm from any other photon candi- 
date or track and have a lateral energy profile consistent 
with a photon shower. We veto photons from a 7r'^{ri) by 
requiring that the invariant mass of the combination with 
any other photon of energy greater than 50 (250) MeV not 
be within the range 115(508) < M^^ < 155(588) MeV/c^. 

The K* is reconstructed from if +, if^, 7r~ and 7r° can- 
didates through the four modes K*° if+7r~, ifgTr^ 
and K*+ K+Tr°, K°tt+ . The K+ and tt" track can- 
didates are required to be well reconstructed in the drift 
chamber and to originate from a vertex consistent with 
the e+e" interaction point (IP). The if^ candidates are 
reconstructed from two oppositely-charged tracks com- 
ing from a common vertex displaced from the IP by at 
least 0.2 cm in the transverse plane and having an in- 
variant mass 489 < M^+^- < 507 MeV/ c^. A track is 
identified as a kaon if it is projected to pass through 
the fiducial volume of the particle identification detector, 
an internally-reflecting ring-imaging Cherenkov detector 
(DIRC) [Q , and the cone of Cherenkov light is consistent 
in time and angle with a kaon of the measured track mo- 
mentum. A charged pion is identified as a track that is 
not a kaon. The tt" candidates are reconstructed from 
pairs of photons, each with energy greater than 30 MeV, 
and are required to have 115 < Mj^ < 150 MeV/ c'^ and 
Et^o > 200 MeV. A mass-constraint fit to the nominal 
tt" mass is used to improve the resolution of its momen- 
tum. The K* reconstruction is completed by requiring 
the invariant mass of the candidate pairs to be within 
100 MeV/c2 of the K*°/K*+ mass. 

The B meson candidates are reconstructed from the 
K* and 7 candidates. The background is predominantly 
from continuum qq production, with the high-energy pho- 
ton originating from initial-state radiation or from tt*^ 
and 77 decays. The background from other B meson de- 
cays is found to be negligible from Monte Carlo simu- 
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FIG. 1: The event shape variable Icos^rj for B° -ft'*°7, 
K*" — > K^n~ Monte Carlo and off-resonance data. 



lation. We exploit event topology differences between 
signal and background to reduce the continuum contri- 
bution. We compute the thrust axis of the event exclud- 
ing the B meson daughter candidates. Figure shows 
the distribution of | cos 9^ \ for signal Monte Carlo events 
and off-resonance data, where 9^ is the angle between 
the high-energy photon candidate and the thrust axis. 
In the center-of-mass frame, BB pairs are produced ap- 
proximately at rest and produce a uniform | cos 9^ \ distri- 
bution. In contrast, qq pairs are produced back-to-back 
in the center-of-mass frame which results in a | cos 9^ \ 
distribution peaking at 1. We require |cos0^| < 0.8. 
We further suppress backgrounds using the angle of the 
B meson candidate's direction with respect to the beam 
axis, 9%, and the helicity angle of the K* decay, 9h- The 
hclicity angle is defined as the angle between either one 
of the K* daughters' momentum vectors computed in the 
rest frame of the K* and the K* momentum vector in 
the parent B meson rest frame. It follows a sin^ 9h dis- 
tribution for the signal and peaks slightly towards ±1 for 
qq background. The B meson candidate's direction also 
follows a sin^ 9"^ for the signal and is approximately flat 
for the qq background. We require | cos 9'^ \ < 0.80 and 
cos 6*^1 < 0-75. 

Since the B mesons are produced via e+e~ T(4S') — > 
BB, the energy of the B meson in the center-of-mass 
frame is the beam energy, E^^^^. This is compared to 
the measured energy of the B meson daughters by defin- 



ing AE* = e*j^,+e; 



- E* 

beam ' 



The distribution of AE* is 



peaked at zero for the signal with a width dominated by 
the resolution of the photon candidates. It is asymmetric 
due to energy leakage from the calorimeter. We require 
-200 < AE* < 100 MeV for the K+n-, K'^Tr+ modes 
and -225 < AE* < 125 MeV for the modes containing 
a 7r°, namely K^tt'^ and KItt^ . The beam-energy sub- 
stituted mass is defined as mEs — y^^^beam^^B ^ where 
p|j is the momentum vector of the B meson candidate 
calculated from the measured momenta of the daughters. 



The mES distribution for the signal is well described by 
an asymmetric resolution function with an approx- 
imately Gaussian core dominated by the resolution of 
the beam energy measurement, and an asymmetric tail 
caused by the energy leakage from the calorimeter for 
the photon candidates. For the modes containing a sin- 
gle photon candidate, namely K~^t:~ and i^gTr^, we can 
remove the tail in toes by rescaling the measured photon 
energy E* by a factor k, determined for each event, so 
that E*i^, +kE*-EI^^^ = 0. The signal for these modes 
is then described by a Gaussian. The background is em- 
pirically described by a threshold function for each 
mode. We select candidates with toes > 5.2GeV/c^. 




FIG. 2: TTiES for the B — > K*'y candidates. The fits used to 
extract the signal yield are described in the text. 



Figure H shows the toes distribution for each of the four 
modes. An unbinned maximum-likelihood technique is 
used to fit the toes distributions for signal j|] and back- 
ground 1 10 contributions. The signal mean and width 
are allowed to vary in the fit for the high statistics K^ti^ 
and if+TT^ modes. The fitted width is slightly larger than 
the predicted Monte Carlo value. In the lower statistics 
modes we fix the width to the Monte Carlo value adjusted 
for the small difference observed in the high statistics 
modes. We fit the on-resonance data with a signal plus 
background shape, and simultaneously the on-resonance 
sideband and off-resonance samples with the same back- 
ground function, using a common fit parameter. The 
off-resonance data sample is required to pass the same 
selection criteria as the on-resonance data sample except 
that we remove the kaon particle identification require- 
ment to gain statistics in the if+Tr" and K'^tt^ modes. 
The on-resonance sideband sample is selected with the 
same criteria as the on-resonance data sample, except 
that we require 150 < AE* < 400 MeV in the X°7r" and 
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iT+Tr" modes, and 100 < l^E* < 500 MeV in the K+n- 
and -ftTgTT"'" modes. The signal yields with statistical er- 
rors from the fit are given in Table ffl. 



TABLE I: The fitted signal yield, efRciency (including B{K*) 
and B{K°)), cross-feed and down-feed from other penguin 
decays and measured branching fraction B{B K'y) for 
each of the decay modes. 



Mode 


Effici 


#Signal 


#Cross 


#Down 


B{B K'-y) 




-cncy 


events 


-feed 


-feed 


ibstat. lb sys. 




% 




events 


events 


Xl0~^ 




14.0 


135. 7± 13.3 


0.4 ± 0.1 


0.6 ± 0.1 


4.24 ± 0.41 ± 0.22 




1.4 


14.8 ± 5.6 


0.4 ± 0.1 


1.0 ± 0.2 


4.10 ± 1.71 ± 0.42 




3.9 


28.1 ± 6.6 


0.7 ± 0.2 


1.2 ± 0.2 


3.01 ± 0.76 ± 0.21 




4.3 


57.6 ± 10.4 


1.2 ± 0.2 


2.6 ± 0.4 


5.52 ± 1.07 ± 0.38 



As a consistency check we plot in Figure ^ the Ai?* 
projection for the K^tt^ mode after requiring 5.27 < 
Toes < 5.29GeV/c^. A comparison of the observed /S.E* 
distribution with Monte Carlo shows good agreement. 
We also plot Mk+-k- in Figure after requiring 5.27 < 
TOES < 5.29GcV/c2, -200 < A^;* < 100 MeV and 0.7 < 
^^K+TT- < l.lGeV/c^. We fit with a relativistic Breit- 
Wigner plus linear background shape and determine that 
the signal is consistent with coming from the K*{892). 




0.1 0.2 0.3 
AE*(GeV) 



FIG. 3: a) The AE* projection for B° K*°-y, K*° 
K^-K~ candidates. The curve is the Monte Carlo expecta- 
tion with a linear background, b) The Mjf+^- projection for 
B" K'^-y, K*° K+TT- candidates with the M^+^- mass 
cut relaxed. The curve is a fit to a relativistic Breit-Wigner 
function with linear background. 

The efficiency for the selection of i? ^ K*j candidates 
is given in Table |. The branching fraction is determined 
from the yield, the efficiency and the total number of 
BB events in the sample. The cross-feed from the other 
B — > K*j modes and the down-feed from B — > Xg^ 
are estimated with Monte Carlo assuming the measured 
branching fractions from the CLEO collaboration (7[ ||] 
for each mode and subtracted from the signal yield. 

The total systematic error is the sum in quadrature 
of the components shown in Tabic The systematic 
uncertainty in the signal yield derives from uncertainties 
in the signal line shape, and cross-feed and down-feed 



contributions. The uncertainty in the signal line shape 
results from the ruES width difference described above. 
To gain statistics in the off-resonance data sample used 
to fit the background function for the K^tt~ and i^+7r° 
modes we relax the kaon identification requirement and 
consequently assign a systematic uncertainty to the as- 
sumption that the background shape is unaffected. The 
error in the assumed branching fractions and final-state 
modeling for B — > Xgj Q gives a systematic error in the 
estimated down-feed from these modes. The tracking effi- 
ciency is computed by identifying tracks in the silicon ver- 
tex detector and observing the fraction that is well recon- 
structed in the drift chamber. We estimate the K^l effi- 
ciency uncertainty by comparing the momenta and flight- 
distance distributions in data and Monte Carlo. The 
kaon identification efficiency in the DIRC is derived from 
a sample of D*+ D°Tr+,D° -> K-Tr+ decays. The 
photon and 7r° efficiencies are measured by comparing the 
ratio of events iV(T=^ h^i:°)/N{T^ ^ h^TT^n^) to the 
previously measured branching ratios [pd]| . The photon 
isolation and 71*^/77 veto efficiency are dependent on the 
event multiplicity and are tested by "embedding" Monte 
Carlo-generated photons into both an exclusively recon- 
structed B meson data sample and a generic B meson 
Monte Carlo sample. The AE* resolution is dominated 
by the photon-energy resolution so that uncertainties in 
the calorimeter energy resolution and overall energy-scale 
cause an uncertainty in the efficiency of the AE* require- 
ment. The photon-energy resolution is measured in data 
using 77*^ and rj meson decays and e+e^ — > e^e~7 events. 
The energy scale uncertainty is estimated by using a sam- 
ple of T] meson decays with symmetric energy photons; 
the deviation in the reconstructed rj mass from the nom- 
inal T] mass provides an estimate of the uncertainty in 
the measured single photon energy. 

The B K*^ samples, except for the K^^tt^ sample, 
are used to search for CP-violating charge asymmetries 
Acp, defined by 



A 



CP 



T{B K*-f) - T{B K*-f) 
r{B -> K*-f) + T{B K*-f) 



The fiavor of the underlying b quark is tagged by the 
charge of the or K*^ in the decay. The on-resonance 
sample for each mode is divided into two CP-conjugate 
samples and the signal yield for each is extracted with 
the same fitting technique as for the branching frac- 
tion measurements. In the fit the background shape 
and normalization, as well as the signal peak and width 
are constrained to be the same for both CP-conjugate 
samples. The measured asymmetries and the asymme- 
try of the background in the sideband regions defined by 
-200 < AE* < 100 MeV, 5.2 < toes < 5.27GeV/c2 are 
given in Table III. 



The systematic uncertainty in the asymmetry is due to 
possible detector effects that cause a different reconstruc- 
tion efficiency for the two CP conjugate decays. This un- 
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TABLE II: The systematic uncertainties in the measurement 
of S(B ^ K'j). 

% Uncertainty in 
B{B K'-y) 













772ES Lill6 sll£Lp6 




7.4 


1.7 


1.9 


Background shape 


1.0 






3.8 


Down-feed modeling 


1.0 


1.5 


1.0 


1.2 


/n^ trackine; efficiency 


2.4 




1.2 


1.3 


_£ 1. Q ijiiii^iC/iiv^ y 




4.5 


4.5 




Kaon identification 


0.7 






1.0 


Photon efficiency 


1.3 


1.3 


1.3 


1.3 


Photon distance cut 


2.0 


2.0 


2.0 


2.0 


7r° efficiency 




2.5 




2.5 


/r] veto 


1.0 


1.0 


1.0 


1.0 


Energy resolution 


2.5 


2.5 


2.5 


2.5 


Energy scale 


1.0 


1.0 


1.0 


1.0 


MC statistics 


1.9 


2.4 


1.5 


2.1 


B counting 


1.6 


1.6 


1.6 


1.6 


Total 


5.3 


10.3 


6.7 


7.0 



The measured Acp in signal and background 



TABLE III: 
samples. 



Mode 


ylcp (signal) 


Acp (background) 




(zbstat. ± sys.) 


(±stat.) 




-0.049 ± 0.094 ± 0.012 


-0.011 ± 0.104 




-0.190 ±0.210 ±0.012 


-0.080 ± 0.080 




0.044 ±0.155 ±0.021 


-0.022 ±0.105 



certainty has been estimated in data with a number of 
known charge-symmetric processes to be less than 2%. 

Finally, we combine the measured branching frac- 
tions for the individual modes using a weighted average, 
B{B+ K*+j) ^ (3.83±0.62(stat.)±0.22(sys.))xl0-5, 
bIb" = (4.23±0.40(stat.)±0.22(sys.))xl0-5. 

The weighted average of the measured CP-violating 
charge asymmetries is Acp{B K*j) — —0.044 ± 
0.076(stat.) ± 0.012(sys.). We constrain -0.170 < 
Acp{B K*-i) < 0.082 at 90% C.L. 
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